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ABSTRACT 

GRB 051022 was undetected to deep limits in early optical observations, but precise astrome- 
try from radio and X-ray showed that it most likely originated in a galaxy at z ^ 0.8. We report 
radio, optical, near infra-red and X-ray observations of GRB 051022. Using the available X-ray 
and radio data, we model the afterglow and calculate the energetics of the afterglow, finding it 
to be an order of magnitude lower than that of the prompt emission. The broad-band modeling 
also allows us to precisely define various other physical parameters and the minimum required 
amount of extinction, to explain the absence of an optical afterglow. Our observations suggest 
a high extinction, at least 2.3 magnitudes in the infrared (J) and at least 5.4 magnitudes in the 
optical (U) in the host-galaxy restframe. Such high extinctions are unusual for GRBs, and likely 
indicate a geometry where our line of sight to the burst passes through a dusty region in the host 
that is not directly co-located with the burst itself. 
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Dark gamma-ray bursts (GRBs) — at the most 
basic level those without optical afterglows — are 
a long-standing issue in GRB observations. Al- 
though in many cases the non-detection of an af- 
terglow at optical wavelengths may simply be due 
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to an insufRciently deep search, or one which takes 
place at late times (e.g. Fynbo et al. 2001), a 
subset of GRBs with bright X-ray afterglows re- 
mains undetected despite prompt and deep optical 
searches (e.g. Groot et al. 1998) and directly im- 
plies suppression of the optical light. 

There arc several plausible explanations for 
this, the most likely being that the burst is at 
high redshift, such that the Ly-alpha break has 
crossed the passband in question, or that there 
is high extinction in the direction of the GRB. 
Examples of both have been found, with a small 
number of GRBs at z > 5 appearing as V and R 
band dropouts (e.g. Jakobsson et al. 2006b; Haislip 
et al. 2006) and some GRB afterglows appearing 
very red at lower redshift, due to effects of extinc- 
tion (e.g. Levan et al. 2006; Rol et al. 2007). 

Identification of GRBs at very high redshifts is 
the key to using them as cosmological probes. The 
proportion of bursts exhibiting high dust extinc- 
tion is also interesting from the point of view of 
estimating the proportion of star formation that 
is dust enshrouded, as well as understanding the 
environments which favor GRB production (Tren- 
tham et al. 2002; Tanvir et al. 2004). 

The detection and follow-up of dark bursts at 
other wavelengths is essential, as it enables 1) the 
modeling of the afterglow, deriving estimates of 
the extinction and energies involved, potentially 
providing information about the direct burst en- 
vironment, 2) pinpointing the burst position in the 
host, to enable late-time high resolution imaging 
and the detection of dust enhanced regions in the 
host, and 3) determination of the properties of the 
GRB host itself, such as the SFR and average host- 
galaxy extinction. 

The High Energy Transient Explorer 2 mission 
(HETE-2; Ricker et al. 2003) detected and located 
an unusually bright gamma-ray burst (Olive et al. 
2005) with its three main instruments, the French 
Gamma Telescope (FREGATE), the Wide field X- 
ray monitor (WXM) and the Soft X-ray Camera, 
(SXC), on October 22, 2005. A 2.5 arcminute lo- 
calization was sent out within minutes, enabling 
prompt follow-up observations (e.g. Torii 2005; 
Schaefer 2005); a target-of-opportunity observa- 
tion was also performed with Swift. Details of the 
HETE-2 observations can be found in Nakagawa 
et al. (2006). 



The Swift observations resulted in the detec- 
tion of a single fading point source inside the SXC 
error region, which was consequently identified as 
the X-ray afterglow of GRB 051022 (Racusin et al. 
2005a). However, optical and near infra-red (nIR) 
observations failed to reveal any afterglow to deep 
limits, while radio and millimeter observations 
with the Very Large Array (VLA), the Wester- 
bork Synthesis Radio Telescope (WSRT) and the 
Plateau de Bure Interferometer detected the ra- 
dio counterpart (Cameron & Frail 2005; Van der 
Horst et al. 2005; Bremer et al. 2005). The posi- 
tion coincides with its likely host galaxy (Berger 
& Wyatt 2005) at a redshift of z = 0.8 (Gal- Yam 
et al. 2005). 

In this paper, we describe our X-ray, optical, 
nIR and radio observations of GRB 051022. The 
outline of the paper is as follows: in Section 2 we 
describe our observations, data reduction and ini- 
tial results. In Section 3, we analyze these results 
and form our afterglow picture, which is discussed 
in Section 4. Our findings are summarized in Sec- 
tion 5. 

In the following, we have used F oc i/~^t~°' in 
our definition of a and /?. We assume a cosmology 
with Ho = 71 kms-iMpc"\ Qm = 0.27 and f^A = 
0.73. All quoted errors in this paper are 1 sigma 
(68%) errors. 

2. Observations and data reduction 
2.1. X-ray observations 

X-ray observations were performed with the 
Swift X-Ray Telescope (XRT) and the Chandra 
X-ray Observatory (CXO). 

The XRT started observing the afterglow of 
GRB 051022 3.46 hours after the HETE-2 trigger, 
for a total effective integration time of 137 ks be- 
tween October 22 and November 6. 

Observations were performed in Photon Count- 
ing (PC) mode, the most sensitive observing 
mode. We reduced the data using the Swift soft- 
ware version 2.6 in the HEAsoft package version 
6.2.0. Data were obtained from the quick-look 
site and processed from level 1 to level 2 FITS 
files using the xrtpipeline tool in its standard con- 
figuration. The first two orbits (until 2.1 x 10^ 
seconds post burst) show pile-up and were there- 
fore extracted with an annular rather than circular 
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region, with an inner radius of 19 and 12" for or- 
bits 1 and 2, respectively, and an outer radius of 
71". Orbits 3-7 (2.4 x 10* - 4.9 x 10* seconds) 
were extracted with a circular region of 71" radius, 
and later orbits were extracted using a 47" radius 
circle instead. The data for the light curve were 
extracted between channels 100 and 1000, corre- 
sponding to 1 and 10 keV, respectively; while the 
commonly used range is 0.3 - 10 keV, the large ab- 
sorption prevents the detection of any data from 
the source below 1 keV. Otherwise, the procedure 
is similar to that described in Evans et al. (2007). 

Observations with the CXO started on Octo- 
ber 25, 2005, 21:14:20, 3.34 days after the HETE 
trigger, for a total integration time of 20 ks (Pa- 
tel et al. 2005). Data were reduced in a standard 
fashion with the CIAO package. 

We performed astrometry by matching X-ray 
sources with an optical ii-band image that was 
astrometrically calibrated to the 2MASS catalog. 
Our CXO position is RA, Dec = 23:56:04.115, 
+19:36:24.04 (J2000), with positional errors of 
0.33" and 0.12" for the Right Ascension and Dec- 
lination, respectively. This puts the afterglow 
within 0.5" of the center of its host galaxy. 

We modeled the XRT spectra with an absorbed 
power law in XSpec (Arnaud 1996), using data 
from the first seven orbits. A good fit (x^/d.o.f . = 
87.2/99) was obtained with a resulting spectral en- 
ergy index of /? = 1.00 ±0.12 and excess absorp- 
tion (at z = 0.8 and for assumed Galactic abun- 
dances) of TVh = (2.82 ± 0.46) x 10^2 cm'^ on 
top of the estimated Galactic absorption at this 
position (A^H = 4.06 x 10^° cm'^. Dickey & Lock- 
man 1990). The CXO data are fully in agreement 
with these values, showing no change in the spec- 
trum over time between 0.3 and 3.3 days after the 
burst. The absorption measured is far less than 
that measured by the HETE team in their prompt 
data, A^H = (8.8l i;^) x 10^^ cm'^ (Nakagawa et al. 
2006). This could indicate a change in absorption 
between the early (prompt) measurements and 
those at the time of the XRT observations. For 
the prompt emission spectriim, however, the val- 
ues found by Konus-Wind (Golenetskii et al. 2005) 
are rather different than those found by HETE-2, 
and may be the result of the lower energy cut-off 
for FREGATE compared to Konus-wind. Alterna- 
tively, the fact that these spectra are an average 
over the whole emission period may also result in 



incorrect model parameters. In the two last cases, 
the A^H in the prompt emission could be as low as 
the XRT value and still produce an equally well 
fit, but with slightly different model parameters. 

For the XRT data, Butler et al. (2005a) and 
Nakagawa et al. (2006) find a value somewhat 
higher than our value (4.9 x 10^^ cm^^ and 5.3 x 
10^2 cm~2 respectively, when scaled by (1 + zf" , 
Gunn & Peterson 1965). This difference could be 
explained by a different count-binning or an up- 
dated XRT calibration used in our modeling. 

The XRT light curve count rates have been con- 
verted to 1-10 keV fluxes using the results from 
our spectral modeling and calculating the ratio of 
the flux and count rate at the logarithmic center of 
the orbits. The 1-10 keV CXO flux was derived 
using the actual spectral fit. 

A broken power law fit to the X-ray light curve 
results in ai = 1.16 ± 0.06, a2 = 2.14 ± 0.17 and 
a break time of 110^23 aroimd 1.27 days. 

The difference between a\ and a2, and the fact 
that the spectral slope does not change across the 
break (the CXO measurement is past the break), 
are highly indicative that the observed break in 
the light curve is a jet break. In Section 3.1, we 
perform full modeling of the afterglow using the 
fireball model, indeed resulting in a jet-break time 
that agrees reasonably well with the break time 
as determined from only the X-rays. We point out 
that our value for is different than that cited in 
Racusin et al. (2005b), largely because their mea- 
surement of was based on a preliminary XRT 
light curve. 

2.2. Optical and near infra-red observa- 
tions 

Observations were obtained in Z and i?-band 
with the William Herschel Telescope (WHT) us- 
ing the Auxiliary Port and the Prime Focus Imag- 
ing Camera, respectively, in r'i' z' with the Gemini 
South telescope using the GMOS instrument, in 
JHKs with the Wide Field Camera on the United 
Kingdom InfraRed Telescope (UKIRT), in BVRI 
with the DFOSC instrument on the Danish 1.54m 
telescope and in J and Ks with the Southern 
Astrophysical Research (SOAR) telescope using 
OSIRIS. The optical data were reduced in a stan- 
dard fashion using the ccdproc package within the 
IRAF software (Tody 1986), whereas the SOAR 
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data were reduced using the cirred package within 
IRAF. The UKIRT data were reduced using the 
standard pipehne reduction for WFCAM. 

Photometric calibration was done using the cal- 
ibration provided by Henden (2005) for Johnson- 
Cousins filters. For the r'i'z' CMOS filters, we 
converted the magnitudes of the calibration stars 
provided by Henden to the Sloan filter system us- 
ing the transformations provided by Jester et al. 
(2005), and verified by the published CMOS zero 
points. The WHT Z-band was calibrated us- 
ing the spectroscopic standard star SP2323-I-157. 
Calibration of the infrared JHK magnitudes waw 
done using the 2MASS catalog (Skrutskie et al. 
2006). 

No variable optical source was found at the po- 
sition of the X-ray and radio afterglow. For the 
early epoch images (< 1 day post burst), we esti- 
mated a limiting magnitude by performing image 
subtraction between this and a later image using 
the ISIS image subtraction package (Alard 2000). 
To this end, artificial low signal-to-noise sources 
were added onto the images, with a Gaussian 
PSF matched in size to the seeing (some artificial 
sources were added on top of existing sources, e.g. 
galaxies, some on the background sky) . We deter- 
mined our upper limit to be the point where we 
could retrieve 50% of the artificial sources in the 
subtracted image. This assumes that the change 
in brightness of any point source on top of the host 
galaxy is sufficient to be seen in such a subtracted 
image. With the difference in time between the 
epochs, this seems a reasonable assumption (for 
example, for a source fading with a shallow power 
law like slope of F oc t^^"'' , the magnitude differ- 
ence between the two WHT Z-band observations 
is « 0.6 magnitudes). 

Photometry of the host galaxy has been per- 
formed using aperture photometry, with an aper- 
ture 1.5 times the seeing for each image, estimated 
from the measured FWHM of the PSF for point 
sources in the images. 

Table 1 shows the log of our optical/nIR obser- 
vations, while Table 2 shows the upper limits for 
any optical/nIR afterglow. 

2.3. Radio observations 

Radio observations were performed with the 
WSRT at 8.4 GHz, 4.9 GHz and 1.4 GHz. We 



used the Multi Frequency Front Ends (Tan 1991) 
in combination with the IVC-f-DZB back end^ in 
continuum mode, with a bandwidth of 8x20 MHz. 
Gain and phase calibrations were performed with 
the calibrators 3C 286 and 3C 48, although at 
one 8.4 GHz measurement 3C 147 was used. Re- 
duction and analysis were performed using the 
MIRIAD software package-^. The observations are 
detailed in Table 3. In our modeling described in 
section 3.1 we have also used the VLA radio de- 
tection at 8.5 GHz from Cameron & Frail (2005). 

3. Analysis 

3.1. Broadband modeling 

We have performed broadband modeling of the 
X-ray and radio measurements, using the meth- 
ods presented in van der Horst et al. (2007). In 
our modeling we assume a purely synchrotron ra- 
diation mechanism. 

The relativistic blastwave causing the after- 
glow accelerates electrons to relativistic veloci- 
ties, which gives rise to a broadband spectrum 
with three characteristic frequencies: the peak fre- 
quency y^, corresponding to the minimum energy 
of the relativistic electrons that arc accelerated 
by the blastwave, the cooling frequency v^-, corre- 
sponding to the electron energy at which electrons 
lose a significant fraction of their energy by radia- 
tion on a timescale that is smaller than the dynam- 
ical timescale, and the self-absorption frequency 
Va., below which synchrotron self-absorption pro- 
duces significant attenuation. The broadband 
spectrum is further characterized by the specific 
peak fiux F,y^max and the slope p of the electron 
energy distribution. 

The dynamics of the relativistic blastwave de- 
termine the temporal behavior of the broadband 
synchrotron spectrum, i.e. the light curves at 
given frequencies. At first the blastwave is ex- 
tremely relativistic, but is decelerated by the sur- 
rounding medium. When the Lorentz factor F of 
the blastwave becomes comparable to 6^^, where 
0-^ is the opening angle of the jet, the jet starts to 
spread sideways. At that time, fj, the temporal 
behavior of the broadband spectrum changes (see 
e.g. Rhoads 1997). 



^See sect. 5.2 at http://www.astron.nl/wsrt/wsrtGuide/node6.html 
^ http: // www. atnf . csiro.au /computing/ software /miriad 
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Table 1 

Overview of optical observations 



Start date 


AT (average) 


exposure time 


filter 


seeing 


telescope & instrument 




(days) 


(seconds) 




(arcsec) 




2005-10-22T23:25:14 


0.4287 


1800 


Z 


0.8 


WHX + API 


2005-10-23T00:22:33 


0.4684 


1620 


J 


1.2 


SOAR + OSIRIS 


2005-10-23X00:56:00 


0.4917 


1620 


Ks 


1.3 


SOAR + OSIRIS 


2005-10-23T00:48:03 


0.5144 


1920 


i' 


0.6 


Gemini South + GMOS 


2005-10-23T01:07:53 


0.5288 


1920 


r' 


0.6 


Gemini South + GMOS 


2005-10-23T01:27:46 


0.5426 


1920 


z' 


0.5 


Gemini South + GMOS 


2005-10-23T06:31:03 


0.7525 


720 


J 


1.4 


UKIRX + WFCAM 


2005-10-23T06:36:39 


0.7526 


360 


H 


1.3 


UKIRX -1- WFCAM 


2005-10-23X06:47:59 


0.7604 


360 


K 


1.3 


UKIRX -t- WFCAM 


2005-10-23T21:15:57 


1.3389 


1200 


Z 


1.0 


WHX + API 


2005-10-24X09:35:10 


1.8467 


720 


K 


0.3 


UKIRX + WFCAM 


2005-10-25X01:34:03 


2.5181 


1602 


Ks 


1.3 


SOAR + OSIRIS 


2005-10-25X02:13:18 


2.5454 


720 


J 


1.2 


SOAR + OSIRIS 


2005-10-25X02:22:02 


2.5698 


1920 


r' 


1.1 


Gemini South + GMOS 


2005-10-25X02:39:59 


2.5792 


1440 


z' 


1.2 


Gemini South + GMOS 


2005-10-26X00:36:58 


3.4785 


1800 


R 


1.4 


WHX-I-PFIP 


zUUo-lU-zO 1 Uz.4o:UD 


o.ODyo 


ouu 


Gunn i 


1 A 
1.4 


Ui\1.04in -|- UrUoL/ 


2005-10-26X03:23:35 


3.5942 


600 


R 


1.9 


DK1.54m + DFOSC 


2005-10-27X01:01:04 


4.4952 


600 


B 


2.3 


DK1.54m + DFOSC 


2005-10-27X02:59:20 


4.5773 


600 


R 


1.6 


DK1.54m + DFOSC 


2005-10-27X02:00:48 


4.5367 


600 


V 


1.8 


DK1.54m + DFOSC 


2005-10-28X02:18:38 


5.5491 


600 


i 


1.4 


DK1.54m + DFOSC 


2005-10-30X02:32:59 


7.5590 


600 


B 


1.8 


DK1.54m + DFOSC 


2005-10-30X04:18:30 


7.6323 


600 


U 


1.8 


DK1.54m + DFOSC 


2005-10-30X01:33:57 


7.5180 


600 


V 


1.4 


DK1.54m + DFOSC 


2005-10-31X03:19:05 


8.5910 


600 


B 


1.0 


DK1.54m + DFOSC 


2005-10-31X01:03:40 


8.4970 


600 


R 


1.0 


DK1.54m + DFOSC 


2005-10-31X02:10:02 


8.5431 


600 


V 


1.0 


DK1.54m + DFOSC 


2005-11-01X01:52:57 


9.5312 


600 


R 


0.9 


DK1.54m + DFOSC 


2005-11-02X02:04:47 


10.539 


600 


V 


1.2 


DK1.54m + DFOSC 


2005-11-03X01:10:34 


11.502 


600 


B 


1.2 


DK1.54m + DFOSC 


2005-11-07X01:25:30 


15.512 


600 


Gimn i 


1.4 


DK1.54m + DFOSC 


2005-11-08X01:40:48 


16.523 


600 


Gunn i 


1.4 


DK1.54m + DFOSC 
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Table 2 
Limiting magnitudes 



filter 


limiting magnitude'^ 


AT (average) 


frequency 


specific flux'' 






days 


Hz 




Ks 


> 20.0 


0.4917 


1.40 • 10" 


< 6.82 


J 


> 20.3 


0.4684 


2.40 • IQi"* 


< 12.3 


Z 


> 22.9 


0.4287 


3.43 • 10" 


< 2.66 


z' 


> 23.5 


0.5426 


3.36 • 10" 


< 1.53 


r' 


> 25.3 


0.5288 


4.76 • 10" 


< 0.305 



'^See text for the definition of the limiting magnitude. 

''Specifc fluxes have been corrected for a Galactic extinction value of Eb-v = 0.04 (Schlegel et al. 1998), 
and converted from magnitudes using the calibration by Tokunaga & Vacca (2005) for the JK^ filters; the 
other filters are on the magnitude AB-system (Oke & Gunn 1983) 



Table 3 

Overview of WSRT radio observations 



Start date 


AT (average) 


integration time 


frequency 


specific fiux 




(days) 


(hours) 


(GHz) ' 


(A*Jy) 


2005-ll-04T18:14:24 


13.37 


4.0 


8.5 


38 ± 132 


2005-ll-08T14:19:41 


17.19 


7.0 


8.5 


28 ± 97 


2005-10-23T15:20:10 


1.19 


5.0 


4.9 


281 ± 32 


2005-10-24T15:17:17 


2.22 


6.2 


4.9 


342 ± 34 


2005-10-25T15:12:58 


3.30 


5.4 


4.9 


143 ± 30 


2005-10-28T18:33:08 


6.40 


8.5 


4.9 


91 ± 28 


2005-10-30T18:00:00 


8.32 


5.8 


4.9 


138 ± 28 


2005-ll-01T18:00:00 


10.38 


8.9 


4.9 


169 ± 28 


2005-ll-04T17:31:12 


13.37 


4.6 


4.9 


70 ± 34 


2005-10-25T15:56:10 


3.33 


5.4 


1.4 


8 ± 78 
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We fit our data to six parameters: Vc. fm-, i^a, 
-Fi/,max, P and tj. From these parameters and the 
redshift of the burst, z = 0.8, we can find the 
physical parameters governing the blastwave and 
its surroundings: the blastwave isotropic equiv- 
alent energy Eiso, the jet opening angle 9^, the 
collimation corrected blastwave energy -Ejct, the 
fractional energy densities behind the relativistic 
shock in electrons and in the magnetic field, £e 
and sb respectively, and the density of the sur- 
rounding medium. The meaning of the latter pa- 
rameter depends on the density profile of the sur- 
rounding medium. For a homogeneous circum- 
burst medium, we simply determine the density 
n. For a massive stellar wind, where the density 
is proportional to with R the distance to the 
GRB explosion center, we obtain the parameter 
j4*, which is the ratio of the mass- loss rate over 
the terminal wind velocity of the GRB progenitor. 

Our modeling results are shown in Table 4, for 
both the homogeneous external medium and the 
stellar wind environment. The light curves for the 
best fit parameters are shown in Figure 1. We 
have performed Monte Carlo simulations with syn- 
thetic data sets in order to derive accuracy esti- 
mates of the best fit parameters, which are also 
given in the table. It is evident from the results 
that our six fit parameters are reasonably well con- 
strained in both cases for the circumburst medium. 
The derived physical parameters are also well con- 
strained, except for and eb- The values we 
find for both the isotropic and the collimation cor- 
rected energy, are similar to those found for other 
bursts; this is also true for p. See e.g. Panaitescu 
& Kumar (2001) and Yost et al. (2003). The jet 
opening angle and the density of the surrounding 
medium are quite small, but both not unprece- 
dented. The jet-break time tj is somewhat smaller 
than estimated in Section 2.1, but both estimates 
have relatively large errors, likely because of the 
lack of (X-ray) data around the jet-break time. 

With the absence of optical light curves, it is 
not possible to discriminate between the two dif- 
ferent circumburst media. This is mainly due 
to the fact that the X-ray band lies above both 
z/m and Vc, in which case the slopes of the light 
curves do not depend on the density profile of 
the circumburst medium (even at 0.15 days, back- 
extrapolating from Table 4 results in its value 
being below the X-ray band). The Xred some- 



what better for the stellar wind case, but the ho- 
mogeneous case cannot be excluded. From the 
X-ray light curve, however, one can conclude that 
the density profile of the medium does not change 
between approximately 0.15 and 12 days after the 
burst. If there were a transition from a stellar 
wind to a homogeneous medium, the X-ray flux 
has to rise or drop significantly, unless the den- 
sities are the fine-tuned at the transition point 
(Pe'er & Wijers 2006). From the fact that the 
medium does not change during the X-ray obser- 
vations, one can draw conclusions on the distance 
of the wind termination shock of the massive star: 
if one assumes that the medium is already homoge- 
neous at « 0.15 days, the wind termination shock 
position is at R„ < 9.8 • 10"'^'^ cm (0.32 pc); if the 
circumburst medium is a stellar wind up to « 12 
days after the burst, i?w ^ 1-1 • 10^^ cm (3.7 pc). 

3.2. The non-detection of the optical af- 
terglow 

It is quickly seen that GRB 051022 falls into 

the category of the so-called "dark bursts". Us- 
ing, for example, the quick criterion proposed by 
Jakobsson et al. (2004), we find Pox < -0.05 
at 12.7 hours after the burst using the Gemini r' 
band observation, well below the proposed limit of 
Pox < 0.5. A more precise criterion would com- 
bine the available spectral and temporal parame- 
ters of the X-ray afterglow, allow all valid combi- 
nations, and from that infer the range of possible 
optical magnitudes from the X-rays (see e.g. Rol 
et al. 2005). This is, in fact implied in our previous 
modeling: the modeled specific fluxes correspond- 
ing to the band and epoch of our optical and nIR 
upper limits are listed in Table 5 (see also Table 
2). 

While the values in this table are given for local 
extinction, not K-corrected to z = 0.8, it is imme- 
diately obvious that our ii'-band observations put 
a stringent constraint on the required extinction 
directly surrounding the burst. 

To estimate the amount of local extinction in 
the host galaxy, we have modeled the nIR to X-ray 
spectrum around 0.5 days after the burst, consid- 
ering 3 different extinction curves: those of the 
Milky Way (MW), the Large Magellanic Cloud 
(LMC) and the Smafl Magellanic Cloud (SMC), 
from Pei (1992), with Ry of 3.08, 3.16 and 2.93, 
respectively. 
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Table 4 

Results of broadband modeling for both a homogeneous external medium and a massive 

STELLAR WIND. THE BEST FIT PARAMETERS ARE SHOWN TOGETHER WITH ACCURACY ESTIMATES FROM 

Monte Carlo simulations with synthetic data sets. The characteristic frequencies of the 
synchrotron spectrum and the specific peak flux are given at t]. 



Parameter 


Homogeneous 


Stellar wind 




(lAB+lil) ■ IQi^ Hz 


(2.84+°-32) . 1017 Hz 


t'm(ij) 


(3.50t2-26) . iQii Hz 


(2.90i? ;3) ■ 10" Hz 




(4.56t^-^^) • 109 Hz 


(2.68t?-^^) • 10^ Hz 


-^i/,max(^j ) 






p 


^■"-"J-0.05 




h 


0.961^-28 days 


1.06+°-^} days 




3.39+1'°^ deg 


2.30tj;°9 deg 


E- 


{5.23t\:ll) ■ 10^2 erg 


(28.2t?J °) • 10^2 erg 




(0.917+|]-^^^) • 10^0 erg 


(2.27t^ ?^) • 10^0 erg 




0.247t-«^ 


0.068lt°;g95i 


£b 


(7.631^^3^^). 10-3 


(8.02t^*i^^) • 10-3 


n 


(1.06t?:^I) • 10-' cm-3 




A 




(2.94t^-?8) . 10-2 


-^red 


1.9 


1.5 



'^Thc parameter is a measure for tire density in the case of a stellar wind environment, being the ratio 
of the mass-loss rate over the terminal wind velocity, and here given in units of 10^^ Solar masses per year 
divided by a wind velocity of 1000 km/s (see van der Horst et al. 2007). 



Table 5 

Upper limits compared to model specific flux calculations. The inferred lower limits on 
the extinction are given in the observers frame. The Eb-v values are given for a 
Galactic extinction curve = 3.08), and are for illustrative purposes; see the 

comments at the end of section 3.2. 



filter 


upper limit 


homogeneous density profile 




stellar wind density profile 






modeled specific flux 


extinction 


Eb-v 


modeled specific flux 


extinction 




(MJy) 


(^Jy) 


(mag.) 




(^Jy) 


(mag.) 


Ks 


< 6.82 


93.1 


2.84 


7.74 


57.2 


2.31 


J 


< 12.3 


117 


2.44 


2.71 


74.1 


1.95 


Z 


< 2.66 


103 


3.97 


2.58 


67.8 


3.52 


r' 


< 0.305 


74.5 


5.97 


2.17 


44.4 


5.41 


z' 


< 1.53 


87.7 


4.40 


2.97 


51.9 


3.83 



Eb- 



6.29 
2.16 
2.29 
1.97 
2.59 
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* 
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•on. 


— ■ . . — 
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0.1 1.0 10.0 

Time (days after burst) 



Fig. 1. — Fit results for a homogeneous circumburst medium (left panel) and a massive stellar wind (right 
panel) . The solid and dash-dotted lines are the best model fits, and the dotted and dashed lines indicate the 
predicted rms scatter due to interstellar scintillation; see the appendix for further details. Also included in 
the figure (and modeling) is the reported VLA 8.5 GHz detection (Cameron & Frail 2005, left-most point in 
the 8.5 GHz subplot). 
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For this, wc used the unabsorbcd XRT flux ob- 
tained from the spectral fit to orbits 3-7 (which 
do not contain piled- up data), and fixed the energy 
spectral slope in the X-rays at /3 = 1 (also from the 
X-ray spectral fit). The optical specific fluxes were 
scaled to the logarithmic mid-observation time of 
the X-ray observations with an assumed a = 1.16 
decline. This estimated optical decay is derived 
from the pre-break X-ray decay value, allowing for 
the cooling break between the two wavelength re- 
gions, and averaging the two possible values for 
ctx — ctopt (-0.25 and 0.25). We can further put 
the most stringent constraint on the broken power 
law spectral shape, by setting the spectral break 
just below the X-rays, at 1.8 x 10^'' Hz, which 
follows from our previous broad-band modeling. 
Our results indicate that, for the aforementioned 
extinction curves, a local extinction of Eb-v ~ 7 
(for all three extinction curves) is necessary to ex- 
plain the /sT-band upper limit. 

We can relate the resulting A^h from our X-ray 
spectral fits to any local Eb-v, using the relations 
found in Predehl & Schmitt (1995), Fitzpatrick 
(1985) and Martin et al. (1989) for N{HI)/Eb-v, 
and adjusting the metallicity in our X-ray ab- 
sorption model accordingly. We obtain Eb-v = 
7.5, 1.54 and 0.84 for a MW, LMC and SMC ex- 
tinction curve respectively, with the MW value 
showing the best agreement with our findings for 
optical extinction (both Nakagawa et al. 2006 and 
Butler et al. 2005b find Eb-v values roughly twice 
as high here, for a MW extinction curve only, since 
their A^h estimate is larger than ours). This, ob- 
viously, depends on the assumption that the MW 
(or otherwise, LMC or SMC) extinction curves are 
valid models to compare with our observed data 
here. Since these data happen to originate from 
just one sight line in a galaxy, this may very well 
not be the case. Further, even if the extinction 
curve is correct, the actual value of Ry may be 
rather different for the host galaxy. Finally, the 
Eb-v - -^H relations show a rather large scatter, 
especially at higher column densities, nor is the 
A^H always derived using X-ray spectroscopy. Our 
above results are therefore approximations, which 
are useful to compare with other (GRB host) stud- 
ies, but should be taken with the necessary cau- 
tion. 



3.3. The host galaxy of GRB 051022 

Using the optical data described above, we fit 
the SED of the host of GRB 051022 using the 
HyperZ program^ developed by Bolzonella et al. 
(2000). The photometry of the host has been per- 
formed using apphot within IRAF, in an aperture 
1.5 times the estimated seeing in the different ex- 
posures. The results are reported in Table 6 (see 
also Ovaldsen et al. 2007). The range of photomet- 
ric magnitudes reported in this paper provides one 
of the most complete broadband optical datasets 
of a GRB host galaxy to date. We fit using the 
eight synthetic galaxy templates provided within 
HyperZ at the redshift of the host, and find that 
the host galaxy is a blue compact galaxy of type 
irregular, with a dominant stellar population age 
of K, 20 Myr, similar to other long GRB hosts 
(Christensen et al. 2005). A moderate amount of 
extinction of Ay w 1 mag is required to fit the 
SED, with an SMC-type extinction curve provid- 
ing a best fit, and the luminosity of the host is 
approximately 1.5 (assuming M^^b = —21); 
these findings are in full agreement with Castro- 
Tirado et al. (2006). The amount of extinction 
in the line of sight towards the GRB required to 
suppress the optical light of the afterglow to the 
observed limits is clearly higher than the Ay value 
found from the host SED: Ay =4.4 magnitudes 
towards the GRB, estimated from blueshifting our 
measured (observer frame) z' band extinction to 
z = 0.8. The host galaxy SED extinction is, how- 
ever, an average value derived from the integrated 
colors of the host. 

The host of GRB 051022 is located in a field 
crowded with galaxies of various Hubble types. 
We perform photometry on several galaxies close 
to the GRB host (within 1 arcminute) to investi- 
gate the possibility that the high star formation 
rate seen in the optical (Castro-Tirado et al. 2006 
report an SFR of » 2OM0yr~^) is induced by 
a recent interaction with one of the neighboring 
galaxies. As formation of high mass stars has also 
been observed to occur in dusty regions in merging 
systems (see e.g. Lin et al. 2007), this could help 
to explain the excess optical extinction towards 
GRB 051022. We performed HyperZ fits to these 
galaxies, and find that none of them is well fit by 
a photometric redshift of z « 0.8. Particularly 



^ See http: / /webast. ast.obs-mip.fr/hyperz 
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Table 6 

Measured host galaxy magnitudes 



filter 


magiiiludo 


magiutudo error 


K 


18.40 


0.04 


Ks 


18.36 


0.09 


H 


19.42 


0.09 


J 


19.92 


0.05 




21.41 


0.05 


z' 


21.30 


0.04 


i' 


21.77 


0.01 


r' 


22.04 


0.01 


R 


21.84 


0.09 


V 


22.30 


0.04 


B 


22.75 


0.02 


U 


> 21.3*^ 





"^AB magnitude 
upper limit 



the two galaxies closest to the GRB host galaxy 
are not compatible with a redshift 0.8, and show 
best fits with photometric redshifts of z « 0.2 - 
0.25. Out of the sample of six galaxies close to the 
GRB host we find that four have best-fit photo- 
metric redshifts in the range 0.20 - 0.25, making 
it unlikely that a possible overdensity of galaxies 
near the host galaxy is due to a cluster or galaxy 
group at the host redshift. 

4. Discussion 

The issue of non-detected ( "dark" ) GRB after- 
glows has received significant interest ever since 
the discovery of the first GRB afterglow, start- 
ing with the non-detection of GRB 970828 to very 
deep limits (Groot ct al. 1998; Odcwahn et al. 
1997). For this particular afterglow, its non- 
detection has been attributed to a dust-lane in 
its host galaxy (Djorgovski et al. 2001). Dust ex- 
tinction as the cause of the non-detection of the 
optical afterglow has been inferred in the case of 
several other GRBs, notably those with a precise 
X-ray or radio position, where one can pinpoint 
the afterglow position on top of its host galaxy 
(e.g. GRB 000210, Piro et al. 2002). 

Optical drop-outs due to high redshift will also 



result in dark bursts, but are harder to confirm, 
since it would require at least one detection in a 
red band, to detect the Lya break. Otherwise, it 
becomes indistinguishable from dust extinction. 

Other explanations of afterglow non-detections 
include the intrinsic faintncss of the afterglow. For 
HETE-2 detected GRBs, this has been inferred 
for e.g. GRB 020819 (Jakobsson et al. 2005). For 
Swift bursts, where rapid and accurate X-ray po- 
sitions are often available, this is a major cause 
of non-detections (Berger et al. 2005), largely at- 
tributed to a higher average redshift. 

In our case here, the host galaxy has been de- 
tected at a relatively modest redshift, which al- 
most automatically points to the dust extinction 
scenario. The radio and X-ray detections even al- 
low us to accurately model the necessary amount 
of extinction between us and the GRB. 

4.1. The burst environment 

The issue of the role of dust extinction in the 

lines of sight towards GRBs is still very much 
an open one. While clear signs of dust depletion 
are seen in several afterglow spectra, the Ay val- 
ues that arc predicted from these depletion mea- 
sures are generally much higher than the observed 
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ones, that can be found from the continuum shape 
(SavagUo & Fall 2004). Recently, selected samples 
of GRB afterglows were homogeneously analyzed 
for X-ray and optical extinction, showing dust to 
gas ratios different from Galactic and Magellanic 
cloud values (Starling et al. 2007; Schady et al. 
2007). Galama & Wijers (2001) and Stratta ct al. 
(2004) had already found dust (optical) to gas (X- 
ray) ratios to be lower than the Galactic value (in 
all cases, however, there is a bias in these samples 
to optically and X-ray detected afterglows). Com- 
parison of neutral hydrogen columns and metal- 
licities of afterglow lines of sight with X-ray ex- 
tinction values (Watson et al. 2007) showed that 
the absorption probed by these two wavelength 
regimes is generally located at different positions 
in the host. In all these cases there may be 
significant biases against bursts with low appar- 
ent magnitudes, preventing optical spectroscopy, 
which are hard to quantify. 

In the case of GRB 051022 there is a significant 
discrepancy between the extinction for the host as 
a whole and that along the line of sight to the 
burst, or at least along our line of sight towards 
the burst. This is perhaps not too surprising if 
one assumes, for example, that the burst occurred 
inside a Giant Molecular Cloud (CMC). Jakobs- 
son et al. (2006a) compared the GRB iV(HI) dis- 
tribution to that of modeled GRBs located inside 
Galactic-like GMCs. They found that the two dis- 
tributions are incompatible, and possibly GRBs 
are more likely to occur inside clouds with a lower 
Af(HI), or alternatively, outside the actual GMC. 
(Note that their study concentrates on bursts with 
z> 2, where the Ly-a absorption is visible in the 
optical wavebands; it is also biased towards opti- 
cally detected afterglows). A GMC could there- 
fore actually be positioned in front of the GRB, 
where the required optical and X-ray extinction 
is easily achieved. This agrees with the findings 
by Prochaska et al. (2007), who analyzed several 
GRB-Damped Lyman Alpha spectra and from ob- 
served depletion levels infer that the gas is not 
located directly near the GRB (e.g. its molecu- 
lar cloud) but further out. The specific case of 
GRB 060418 confirmed this through time-resolved 
high resolution spectroscopy, showing that the ob- 
served metal lines originate past 1.7 kpc from the 
burst itself (Vrecswijk et al. 2007). In fact, X- 
ray radiation from the burst could easily destroy 



grains out to 100 pc (Waxman & Draine 2000; 
Fruchter et al. 2001; Draine & Hao 2002) and 
permit the afterglow radiation to penetrate the 
surrounding molecular cloud. Dust extinction is 
therefore likely to occur further out, perhaps to 
several kiloparsecs. 

It is interesting to find a non-SMC type of ex- 
tinction curve from the combination of X-ray and 
optical absorption (though not completely ruled 
out): in most cases modeled, an SMC extinction 
curve fits the optical-X-ray spectra best (Starling 
et al. 2007; Schady et al. 2007), presumably at- 
tributable to the absence of the 2175 A feature 
(Savage & Mathis 1979) and the low dust to gas 
ratio. Our findings indicate that the extinction 
along the line of sight to the GRB will generally 
be different than one of the three assumed extinc- 
tion curves. Local small scale density variations in 
clouds, such as found by from infrared studies in 
the Taurus region and from simulations (Padoan 
et al. 2006), could cause this fairly easily. 

4.2. Energetics 

Our modeling provides us with a detailed set 
of parameters of the afterglow energetics, includ- 
ing -Ejct! the energy of the afterglow. For the 
prompt emission energy, we use the data from 
the Konus-Wind measurements (Golenetskii et al. 
2005). We calculate a prompt isotropic energy of 
4.39i^;i^ X 10^3 erg in the 20 keV - 20 MeV ob- 
server frame, and, by applying a K-correction (as 
in e.g. Bloom et al. 2001), £;p,iso = 10.41:^;^ x 10^^ 
erg in the 1-10^ keV rest frame. The coUimation 
corrected energy depends on the assumed density 
profile of the surrounding medium: for a homoge- 
neous medium, we obtain -Ep,jet = 18-2 x 10^° erg, 
and for a wind-like medium, E'pjet = 8.38 x 10^° 
erg. With E'poak = 9 18^59 keV in the burst rest 
frame, we find that the i^peak - -E^p,jet relation 
(Ghirlanda et al. 2004) somewhat underestimates 
the -Bpcak when calculated from -Ep.jct: i'pcak ~ 
740 keV for a homogeneous medium, and ~ 430 
keV for a wind medium (the difference between 
our chosen cosmology and that used by Ghirlanda 
et al. 2004 amounts to only a 0.3% difference in 
Eiso)- These estimates, however, come with a few 
caveats: 1) the Speak from the Konus-Wind data 
is calculated using an exponential cut-off model, 
not the Band function (Band et al. 1993). Since 
the Band function includes the case of an expo- 



12 



ncntial cut-ofF model (with /? = —oo, this should, 
however, pose no problem in estimating the actual 
-Epeak), 2) our break time, and therefore the jet- 
opening angle, are calculated from the full model- 
ing of the afterglow, which effectively means de- 
rived from the available X-ray and radio data. 
Further, the original Ghirlanda relation was de- 
rived using optical break times. Recent efforts 
show that estimating jet-break times from X-ray 
light curves may not lead to the same results (e.g. 
Panaitescu et al. 2006), and 3) the relatively large 
error on the jet opening angle estimate allows for a 
relatively large range in coUimation corrected en- 
ergies. We have simply used here our best value, 
but an Epeak value of 1498 keV derived from Ejet 
can still be accommodated within our errors. (We 
note that, with a different ii'peak estimate and an 
incorrect value for the jet-break time, Nakagawa 
et al. 2006 still found their results to lie on the 
Ghirlanda relation). The break time problem can 
be avoided by looking only at the Epe^k - ^^p.iso re- 
lation (Amati et al. 2002; Amati 2006). From this, 
we estimate ii'peak ~ 924 keV, nicely in agreement 
with the value found directly from the spectra fit. 

Comparing the prompt emission energy (Ep^^et) 
and afterglow blast wave kinetic energy (i5jet)j we 
find their ratio to be i?pjet/£'jct = 3.7 in the case 
of a wind-like circumburst medium, while for a 
homogeneous medium, iJpjetZ-Ejet = 20. These 
ratios are similar to those found for other bursts 
(e.g. Berger et al. 2003, Figure 3). 

GRB 051022 is also one of the brightest bursts 
observed by HETE, with a prompt 30-400 keV 
fluence oi S = 1.31 x 10~^ erg cm~^ (Nakagawa 
et al. 2006). In fact, compared to the sample of 
35 FREGATE bursts analyzed by Barraud et al. 
(2003), GRB 051022 has the largest fluence, even 
topping the relatively c;lose-by GRB 030329 (Van- 
derspek et al. 2004, S = 1.2 x IQ-* erg cm^^; 
note that for GRB 051022, its redshift is close to 
the median redshift of HETE-2 detected GRBs 
and therefore distance effects will play a very mi- 
nor role). Rol et al. (2005) noted this potential 
correlation of fluence with the non-detection of a 
GRB afterglow for the small subset of genuinely 
dark bursts in their sample: the truly dark bursts 
all have a much higher than average fluence (al- 
though this is for a relatively small sample only). 
Potentially, this could point to an external origin 
for the prompt emission, instead of being due to 



internal shocks: a large amount of dust may re- 
sult in more matter that will radiate, while at the 
same time the radiation will be suppressed at UV 
and optical wavelengths. This would indicate an 
origin of the extinction quite close to the burst in- 
stead, in contrast to previous findings for other 
bursts, as discussed in Section 4.1. These lat- 
ter bursts, however, were all optically selected to 
obtain spectroscopy, and may therefore show dif- 
ferent surroimdings than GRB 051022. Unfortu- 
nately, with the small sample size of genuine dark 
bursts a firm conclusion on this correlation is not 
possible, but remains something to watch for in 
future dark bursts. 

5. Conclusions 

GRB 051022 is a prototypical dark burst, with 
the local extinction exceeding 2.3 magnitudes in 
J and 5.4 magnitudes in U, in the host-galaxy 
restframe, with the exact limits depending on the 
circumburst density profile. The extinction curve 
derived from an X-ray - optical spectral fit points 
towards a Galactic type of extinction curve, al- 
though it is likely that this is more or less a co- 
incidence: the host galaxy itself is best modeled 
with an SMC-like extinction curve, with a mod- 
est amount of extinction, Ay « 1 mag. The 
large optical absorption towards the afterglow of 
GRB 051022 is therefore probably the effect of an 
unfortunate position in the host where the line of 
sight crosses dense regions within the host. 

The X-ray and radio afterglow data allow for 
a full solution of the blastwave model, although 
we unfortunately cannot distinguish between the 
density profile (homogeneous or wind-like) of the 
circumburst medium. We estimate a collimation- 
corrected energy in the afterglow emission of 0.92 - 
2.3 X 10^° erg, while the energy in prompt emission 
(1 - 10^ keV rest frame) is 8.4 - 18 xlO""^" erg. 
Aside from the large optical extinction, the af- 
terglow otherwise appears as an average after- 
glow, with no outstanding properties. The poten- 
tially interesting point here is that the 30-400 keV 
fluence of the prompt emission is one of the largest 
ever detected in the HETE-2 sample. 

In the era of Swift GRBs, dust-extincted bursts 
can actually be found in optical/nIR thanks to the 
rapid availability of precise positions: examples 
are found where the burst is relatively bright early 
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on at optical/nIR wavelengths, while the afterglow Center (HEASARC), provided by NASA's Goddard 

proper (post few hours) often can go undetected Space Flight Center, 

(e.g. Gates et al. 2006; Perley et al. 2007). This 
allows targeted follow-up of such dark bursts, i.e. 
determining the host galaxy (and the bursts pre- 
cise position therein) and a redshift measurement. 
In our case, a precise CXO and radio position pin- 
pointed the host galaxy, but such data may not 
always be available. High resolution late-time ob- 
servations of the host, at the location of the GRB, 
may then reveal whether the burst indeed occurred 
inside a dense host region. 
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A. Interstellar scintillation in the radio modeling 

The 4.9 GHz measurements show scatter around the best fit hght ciirvc, which can be accounted for by 
interstellar scintillation (ISS). In Figure 1 we have indicated the predicted rms scatter due to ISS. We have 
calculated the scattering measure from the Cordes & Lazio (2002) model for the Galactic distribution of 
free electrons: SM = 2.04 ■ 10^^ kpc /m^^"/''. The radio specific fiux will be modulated when the source 
size is close to one of the three characteristic angular scales, i.e. for weak, refractive or diffractive ISS. From 
Walker (1998), we calculate the transition frequency between weak and strong ISS, vq = 9.12 GHz, and the 
angular size of the first Frcsnel zone, 9y„ = 0.994 /ias. Our measurements were all performed at frequencies 
below uq, i.e. in the strong ISS regime, which means that only refractive and diffractive ISS modulate the 
specific flux significantly. We calculate the evolution of the source size in the extreme relativistic phase 
{Of, = R/r) and after the jet-break {9s = ROj), and compare this source size with the diffractive angular scale 
Od = 6-pa{volv)~^^^ = 0.0701 • t'QH^ /uas and the refractive angular scale 9^ = 9Fg{uo/v)^^/^ = 128 • i^qh^^ 
to calculate the modulation index mp. In the case of diffractive ISS the modulation index is 1, and in the case 

of refractive ISS nip = (z/o/i')~^^^^° = 0.286 ■ fQHz°- Because of the expansion of the blastwave the angular 

source size exceeds one of the characteristic angular scales at some point in time. Then the modulation will 
begin to quench as mp{9d/9s) in the case of diffractive ISS, and as 'mp{0^/9sf in the case of refractive ISS. 
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